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ABSTRACT 

Pre-main-sequence (PMS) binaries and multiples enable critical tests of stellar models if masses, metallici- 
ties, and luminosities of the component stars are known. We have analyzed high-resolution, high signal-to-noise 
echelle spectra of the quadruple-star system HD 98800 and using spectrum synthesis, computed fits to the com- 
posite spectrum for a full range of plausible stellar parameters for the components. We consistently find that 
sub-solar metallicity yields fits with lower \ 2 values, with an overall best-fit of [M/H] = —0.20 ± 0.10. This 
metallicity appears to be consistent with PMS evolutionary tracks for the measured masses and luminosities of 
the components of HD 98800 but additional constraints on the system and modelling are needed. 
Subject headings: stars: abundances — stars: individual (HD 98800) — stars: pre-main sequence 



1. THE HD 98800 SYSTEM 

HD 98800 is a multiple system of pre-main-sequence 
(PMS) stars that is particularly well-suited to testing evolu- 
tionary models. It has two visual com ponents dlnnesll 19091) . 
separated by 0.848 arc sec in 2001 dPrato et all l2001h and 
estima ted to have an orbital period of about 300 to 430 
years dTokovininll 19991) . In addition, the brighter (southern) 
A' component is a single-lined spectroscopic binary with a 
period of 262 days, while the fainter 'B' component is a 
double-lined spectroscopic binary w ith a period of 315 days 
dTorres et alj[l99 5; Boden et al. 2005). Thus there are at least 
four stars in the system, designated as Aa, Ab, Ba, and Bb. Of 
these, the light from only three of the stars, Aa, Ba, and Bb, 
is vis ible in a combined optical spectrum (Soderblom et al. 
[19981). 

iPrato et all ((2(301) fitted the estimated luminosity and ef- 
fective temperature of Aa with pre-main sequence evolution- 
ary models and de termined its mass to be 1.1 ± O.1M . 
Bod en et all (120051) measured the orbit of the B-subsystem 
and determined masses of Ba and Bb as 0.699 ± 0.064 and 
0.582 ±0.051 Mq respectively, by using Keck Interferometer 
observations combined with Hubble Sp ace Telescop e (HST ) 
astrometry and radial velocities from [Torres et"al] d!995|). 
Hippa rcos measured a distance o f 46.7 ± 6.2 pc dESAl 
1 1997b . Using their orbital solution, iBoden et all d2005l) de- 
rived an independent distance estimate, placing the system at 
42.2 ± 4.7 p c . A ne w reduction of the Hipparcos data by 
Ivan Leeuwenl d2007bl) gives a parallax of 22.28 ± 2.3 arc- 
sec (see Ivan Leeuwenll2007al) . corresponding to a distance of 
44.9 ± 4.6 pc. Based on the lit hium abundances of the compo- 
nents, [Sodertiomitin d!998l) estimated the age of HD 98800 
to be 10 ± 3 Myr. In addition, HP 98800 is a member of the 
TW Hya association (Kastne r" et al.ll 19971) . 

With well-determined masses and parallax, HD 98800 pro- 
vides a unique opportunity for testing PMS evolutionary mod- 
els. Not only are the masses and luminosities constrained for 
HD 98800, but having several stars in a single system helps 
to distinguish between stellar evolutionary models. This is 
because we assume that the stars in this system are of the 
same age and composition and hence should lie near a sin- 
gle isochrone in a Hertzsprung— Russell (HR) diagram. The 
available model isochrones for low mass stars differ from one 



another significantly, primarily because of the different ways 
the models treat convection. Until now it has not been possi- 
ble to distinguish between the models with any certainty be- 
cause the metallicity of the system was unknown. Further, 
based on the locations of the Ba and Bb components in an 
HR diagram. IBoden et all d2005l) sh owe d that the PMS evolu- 
tionary models of lSiess et alJ d2000h and lBaraffe etalJ (1998) 
yielded masses inconsistent with their dynamical masses, un- 
less HD 98800 has significantly subsolar metallicity of about 
-0.3 to -0.5 dex. 

One might be tempted to adopt for the metallicity of HD 
98800, a value d erived from ano ther member of the TW 
Hya association. lYang et all d2005l) estimated the metallic- 
ity of TW Hya itself, but the result ([M/H] = -0.10 ± 
0.12 ) was imprecise since TW Hya is an accreting T Tauri 
star dMuzerolle et alj|2000t lAlencar & Batalhall2002l) , which 
added uncertainty to the spectrum analysis. Other members 
of the association are much cooler and the available molec- 
ular line data, required for low temperature atmospheres, are 
inadequate for precise abundance analyses. Aside from multi- 
plicity, HD 98800 is a useful system for an abundance analysis 
because its components are warm enough to have spectra sim- 
ilar to the Sun's (i.e., they do not have significant molecular 
features). 

In this paper we put observational constraints on the metal- 
licity of the HD 98800 system from analysis of high resolution 
and high signal-to-noise spectra. Our analysis is complicated 
by that which makes HD 98800 so interesting: the multiplic- 
ity of the system. Components A and B are not far enough 
apart to observe separately from the ground without the aid 
of adaptive optics, and so the spectrum observed consists of 
three stellar spectra added together. The composite nature of 
the spectrum means that our ability to determine the abun- 
dances is lessened because continuum emission from all the 
companions dilutes the depths of the spectrum features. We 
circumvent this by fitting synthetic spectrum models for a va- 
riety of temperatures and radii and find that the best fit models 
for the system indicate subsolar metallicity. 

2. OBSERVATIONS 

We obtained spectra of the HD 98800 star system on three 
nights in 2006 December (Table [TJ with the H amilton echelle 
spectrograph at Lick Observatory IVoglfT 987). when the spec- 
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FIG . 1 . — Observed spectrum (histogram) of HD 98800 in the vicinity of the 
isolated Mn I 6021.793 A line. The three components in the observed spec- 
tram correspond to contributions (dashed curves) from Ba (blue-shifted), Aa 
(central), and Bb (red-shifted). The sum of the components (smooth curve) 
matches the observations well. 



TABLE 1 
Observation Log 



Night 


HJD 


Exposure 
(s) 




H.A. 




S/N 


I 


2454075.062 


1500 


-01 


: 12 : 


47.7 


95-120 


II 


2454076.060 


1500 


-01 


: 12 : 


43.1 


120-160 


III 


2454077.067 


1500 


-00 


: 58 : 


12.2 


110-150 



tral lines of the Aa, Ba and Bb components were visible and 
well separated in radial velocity (see Figure [T]). The data 
were reduced using standard IRAF procedures in the echelle 
package. The resolving power of the instrument was about 
40,000. The spectra have a signal-to-noise ratio per pixel of 
over 100 on each night of observation and cover wavelengths 
from 3920 A to 9920 A. 

At the time of observation the A and B components were 
separated by 0.848 arcsec at position angle 4.4 and the sys- 
tem was observed at hour angles near 1 hour east, meaning 
that A and B were separated by 0.3 arcsec along the disper- 
sion direction at the slit. A very wide slit (1.2 arcsec) was 
used to facilitate the recording of the full light from both the 
visual components. Despite this, we found that the relative 
line depths of Aa to Ba+Bb were not constant from one night 
to the next, which we interpret as a variation in the placement 
of A and B in the slit from night to night. This effect adds 
uncertainty to the Aa/(Ba + Bb) flux ratio, as we note below. 
However, it should not affect the Ba/Bb flux ratio since the 
Ba/Bb separation is negligible compared to the slit width and 
the B subsystem is not resolved. 

3. SPECTRUM ANALYSIS 
3.1. Synthetic spectra 

We fitted the observed spectra by computing separate syn- 
thetic spectra for each of the thr ee stars using Spect r oscop y 
Made Easy (SME), described in IValenti & Piskunovl d 19961) . 
SME is a set of IDL routines coupled to a C++ library that 
contains modules for computing opacities and elemental ion- 
ization fractions as well as radiative transfer. The radia- 
tive transfer code assumes local thermodynamic equilibrium 
throughout the line-forming depths of the stellar atmo sphere. 

Model atmospheres for th e spectrum syn thesis dKuruczl 
1 1992b in SME are taken from iKuruczl d 1993b . who provides 



TABLE 2 

Input parameters for radial velocity estimates 



Component 


log 9 


v sin i 


T c g 


R/Rq 








(km s" 1 ) 


(K) 






Aa 


4.25 


5.0 


4500 


1.75 


0.69 


Ba 


4.21 


3.0 


4050 


1.09 


0.24 


Bb 


4.34 


0.0 


3550 


0.63 


0.07 



opacities on a grid of various stellar parameters. The method 
of interpolating spectra (or model atmospheres, as SME does) 
to find the best fit requires knowledge of the star's effec- 
tive temperature (T c g), surface gravity (logg), metallicity 
([M/H]), microturbulence (£), radial-tangential macroturbu- 
lence (S), projected rotational velocity (v sin i) and radial 
velocity (w r ad)- Here [M/H] is a scale factor applied to the 
solar elemental abundances to generate stellar abundances. 

Throughout our analysis, we used the log q value s for the 
Ba and Bb components from iBoden et alj (120051) (hence- 
forth B05), empirical y sin i values for Aa and Ba from 
Soderblom et alj d 1996b and nominal values of 4.25 for the 
gravity of Aa and 0.0 for the v sin i of Bb (Table |2). Given 
our resolution and signal-to-noise, we cannot distinguish 
v sin i for Bb from zero. W e took £ to be 0.85 km s _1 
from IValenti & Fischer! 42005) (henceforth VF05) and used 
VF 05's relation between 5 and T e g. While previou s stud- 
ies dPadgettlll996T;lJames et all 2006; Sa ntos etall2008h have 
found microturbulent velocities in T Tauri stars in the 2-4 
km s^ 1 range, HD98800 is on ly slightly above the main se- 
quence ( Soderblom et al. 1996) and is not nearly as active as 
most classical T Tauri stars. A higher value for £ than 0.85 
would force our derived metallicity to be even lower than we 
find it, which as we will show below, would only exacerbate 
the already low [M/H] that we find for this system. 

Using the above parameters, we synthesized spectra for 
each of the three stars in the HD 98800 spectrum, then shifted 
each synthesized spectrum by the observed radial velocities 
(Section [3~3l of the stellar components and co-added the three 
shif ted s ynthetic spectra in proportion to the flux ratios (Sec- 
tion [3T2J of the stars to obtain the final modeled spectrum. 



3.2. Computed flux ratios 

Flux ratios for the three stars, which enter into the above 
step, were computed by taking the Planck function, B(X, T) 
at the center of our wavelength range at 6115 A and using 
initial estimates of the radii, Ri of the stars from B05. We 
have for the fractional flux contribution from any one star in a 
spectroscopic multiple system of N stars 



R?B(\,Tj) 



(1) 



with N = 3in our case. The combined spectrum modeled as 
seen from the Earth is then given by 



,v radj ,(v sin 



(2) 



where fj is the intrinsic spectrum of star j (normalized to 
unity) at the rest wavelength Ao and is a function of the 
star's radial velocity (v ra d)j and projected rotation velocity 
(v sin i)j. An example of models for the individual stars 
weighted by their flux ratios and appropriately shifted in ra- 
dial velocity, and the combined spectrum are shown along 
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TABLE 3 

Radial velocity estimates (km s _1 ) 



Component Night I Night II Night III 

o5s. ZP obs! oBs! 3^ 



a Diff erence (model — observed) between model radial velocities ( Torres et al. 
1995) and measured values 

with the corresponding segment of the observations in Figure 
ffl 

3.3. Determination of radial velocities and an improved 
ephemeris 

The radial velocities we measure are necessary input to the 
code as part of the fitting process. We first determined the ra- 
dial velocities of the three stellar components by taking nom- 
inal values of the temperatures and radii for the stars and fol- 
lowed the synthesis procedure outlined above. Table[2]shows 
the parameters of the three stars used to generate spectra for 
radial velocity fitting. These values are close to the final val- 
ues obtained after minimizing the \ 2 statistic between the ob- 
served and modeled spectra (Section l331 >. 

These radial velocity estimates were determined from vi- 
sual inspection, which was adequate to match the models to 
the observed spectrum to determine metallicity. We did not 
base our radial velocity estimates on the \ 2 values for the fits 
since the final mask (Section [3~4l for the combined spectrum 
changes with the radial velocity parameters for the three stars, 
introducing a bias in the behavior of \ 2 ■ The X 2 values be- 
tween different radial velocities are not comparable and hence 
not good estimators of when the radial velocities have been 
best determined. However, based on visual examination, dif- 
ferences in u rac j of 0.2 km s _1 for Aa, 0.5 km s _1 for Ba and 
1.5 km s _1 for Bb were easily discernible and we take these 
as conservative estimates of our measurement precision. 

Table [3] lists the measured radial velocities, along with a 
compa rison to the modeled orbital curve from iTorres et all 
(1995) (henceforth T95), in the sense of model minus ob- 
served. Our measured radial velocities differ significantly 
from the predicted values (Figure As noted above (Sec- 
tion [2), we used a wide slit in order to minimize slit losses. 
This allowed the placement of the stars in the slit to differ 
between nights, causing the zero-point velocity to change. A 
shift of the order of 50% of the slit width would correspond 
to movement of the spectrum by 1 pixel, or 2.5 km s _1 at the 
detector. Further, the spectrum on the CCD shifts vertically, 
depending on the zenith angle of the object to which the tele- 
scope is pointed. This introduces errors into the wavelength 
calibration, since the fixed Th-Ar lamp utilized for this pur- 
pose would form spectral features on the detector at the same 
pixels, independent of the zenith angle. This calibration error 
may be significant for HD98800 (~ 1 to 2 km s _1 ), since 
it was observed at low elevation angles, where the effect is 
greatest. 

However, offsets in the wavelength calibration do not affect 
the precision with which the relative velocities of the compo- 
nents can be measured, since the light from each component 
of the multiple system is displaced by the same amount. Also, 
the position angle of the A-B system varied with respect to the 
slit between nights by 3.6 degrees, and this could, in princi- 



FlG. 2. — Observed minus predicted radial velocities of each component 
(Aa, Ba, Bb) of HD 98800 and each observing night (1,11,111) assuming the 
nominal ephemeris 
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pie, give rise to variations in the relative separation between 
the lines of the A and B systems. But this effect is negligible, 
leading to a relative shift of ~ 0.2 km s _1 . 

We also note that the predicted and observed radial velocity 
separations between Ba and Bb differ systematically by about 
6 km s _1 for all three nights, which is significantly larger than 
expected from our measurement precision (Figure |2). It is 
possible that this is an artifact due to differences in the ve- 
locity scales between this work and the earlier radial veloc- 
ity measurements of T95, on which the predictions are based. 
However, this would imply an unlikely systematic underesti- 
mate of all our radial velocities by about 10% relative to T95. 
As an alternative explanation, we hypothesized that the radial 
velocity predictions were offset from the true orbit due to in- 
accuracies in the orbital period of HD 98800 B. 

To test this, we determined the date when the radial velocity 
separation between Ba and Bb predicted by T95 was closest 
to the observed separation. We found that if we used the pre- 
dictions of T95 for a date 4.5 days later than the date of our 
observation, then the systematic effect vanished and the rms 
separation was reduced to 0.9 kms -1 (Figure [3}. This im- 
plies an orbital period of 314.89 ± 0.03 days, which is 0.26 
days shorter than estimated by T95 and also closer to the value 
of 314.327 ± 0.028 determined by B05. However, this esti- 
mated orbital period is still significantly discrepant from B05. 
Continual measurements of radial velocities covering a wider 
range of epochs are needed to adequately resolve this issue. 

3.4. Matching observed and modeled spectra 

VF05 used SME to determine the Na, Si, Ti, Fe and Ni 
abundances in 1,040 nearby F, G, and K stars. They used 
spectra spanning wavelengths from 5160 A to 5190 A and 
6000 Ato 6180 A, from the Lick, Keck and AAT planet search 
programs. As a prerequisite to their analysis, they tuned the 
atomic line data to match the observed spectrum of the Sun 
with synthesized solar spectra. We used the same wavelength 
interval as VF05, except that we avoided the Mg b-band (5 160 
- 5190 A) due to difficulties with continuum normalization 
in that region. We note that in VF05 this band was used to 
constrain the stellar gravity, a parameter that we do not fit 
for in our analysis. Our wavelength range of 6000 - 6180 A 
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8.80 
-20.60 
36.42 



0.42 
-1.95 
3.24 



7.38 
-20.82 
33.58 



1.89 
-1.35 
5.61 



7.42 
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TABLE 4 
Best fit results 



Night 






^Bb 


^Ba/^Aa 


^Bb/^Aa 


[M/H] 


CAa 


CBa 


CBb 


x 2 




(K) 


(K) 


(K) 
















I 


4535 


4200 


3500 


0.6 


0.4 


-0.21 


0.686 


0.229 


0.085 


9.2 


II 


4535 


4200 


3500 


0.7 


0.5 


-0.18 


0.607 


0.276 


0.117 


16.0 


III 


4515 
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3500 


0.7 


0.5 


-0.21 


0.606 


0.276 


0.118 


15.4 



FIG. 3. — Same as Figurefl] except by using predicted radial velocities 
from the ephemeris 4.5 days after the actual observing dates 
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includes the rest of the range used by VF05. This enabled us 
to adopt line parameters from VF05 for computing our model 
spectra. 

VF05 also found that some spectral features were unsuit- 
able for analysis, certain Ca lines in particular being impossi- 
ble to match over a broad range of temperatures with their 
models. They thus generated wavelength masks that were 
used to define the pixels used to compute the x 2 goodness 
of fit for their modeled spectra. For fitting the HD 98800 
spectrum and determining the x 2 statistic, we modified these 
masks to encompass the lines from the three stars with their 
different radial velocities. This was done by taking the mask 
for a single star, shifting it in radial velocity for each star of 
the 98800 spectrum and then combining the three masks to- 
gether to generate a final mask that included all individual 
'good' pixels from the three masks. This procedure increases 
the number of good pixels, but we still avoided the fitting of 
most 'bad' lines, such as the Fe I line at 6065.482 A, the Fe I 
lines near 6137 A, and the Ca I 6122.2 line, which has broad 
damping wings. Full details of the lines that were fitted and 
the relevant masks in the lab frame (for a single star) can be 
found in VF05. 

We validated our procedure of spectrum fitting using ob- 
servations of 11 nearby G dwarfs with T e g between 5400 
and 6000 K that were modeled earlier by VF05 using SME. 
We observed these stars contemporaneously with our obser- 
vations of HD 98800 using the same instrument and set up 
(the Lick Hamilton spectrograph). Using SME, we found that 
we could match the metallicities found in VF05 within 0.04 
dex, which places our resu lts on the same abundan ce scale as 
VF05. Details are given in lSoderblom et al.1 d2009l) . 



3.5. Determination of abundances 

Having estimated the radial velocities for a given night, we 
proceeded to fit for the metallicities and temperatures of the 
three components. To do this, we computed and compared the 
X 2 goodness of fit between the modeled spectrum and the ob- 
served spectrum for different parameters of the three stars. We 
varied the temperatures and radii (to calculate the flux ratios: 
see equation [TJ of the three components as well as the overall 
metallicity of the system over a grid of values, generating one 
X 2 value for each possible parameter combination. As an ini- 
tial grid, we took five temperature steps of 150 K starting at 
4200 K for Aa, 3750 K for Ba and 3550 K for Bb. Five steps 
were also taken for each of the stellar radii centered on the 
values from B05 with steps equal to the quoted uncertainties 
(i? Aa : 1-27, 1.43, 1.59, 1.75, 1.91 i? Q ; i? Ba : 0.67, 0.81, 0.95, 
1.09, 1.23 R Q ; i? B b: 0.52, 0.63, 0.74, 0.85, 0.96 R Q ). Our se- 
lected [M/H] points were +0.3,0.0, -0.1, -0.2, -0.3, -0.4 
and —0.6 dex. Thus our range of parameter values included 
the temperatures, radii, and metallicities determined by B05. 
We then looked for the model with the minimum x 2 an d 
repeated this process for spectra taken on each of the three 
nights. 

Based on the results of this preliminary analysis with a 
coarse grid, we adopted a finer grid centered on the best-fit 
T c ff and [M/H] values from the coarse-grid analysis. The fine 
grid had five temperature steps, with T Aa ranging from 4495 
to 4575 K (20 K steps), T Ba from 4050 to 4350 K (75 K steps) 
and T Bb from 3500 to 3800 K (75 K steps). We also fixed the 
radius of Aa to unity and normalized the other radii with re- 
spect to i?A a , taking both i?B a /-RA a and Rsb/RAa from 0.0 
to 1.0, thereby obtaining flux ratios for Aa:Ba:Bb from 1:0:0 
to approximately 0.33:0.33:0.33. From the coarse grids, it 
was clear that an [M/H] value near —0.2 dex produced the 
minimum x 2 values, and so we took [M/H] from —0.27 to 
—0.12 dex, spaced at 0.03 dex. The parameter combinations 
giving the minimum x 2 f° r spectra taken on each of the three 
nights with this fine grid are shown in Table [4] 

The Teh values favoring the best fit are at the lowest tem- 
perature in our grid of model atmospheres (3500 K). Our es- 
timate for Teb is also lower than the value (4000 K) in B05. 
We note that the few lines we had to work with for Bb only 
provided weak constraints on its temperature. Therefore we 
adopt 3650 K, which is the next available grid point in tem- 
perature for Bb from the optimal-x 2 solution at 3500 K, as 
an estimate of its temperature. Further, we find that the best 
fit values of T\ a for a given [M/H] are correlated with [M/H] 
itself. This can be understood by the fact that the Aa com- 
ponent has a dominant contribution to the spectrum and that 
increasing metallicity increases the depths of spectral lines, 
which is offset in these fits by a rising best fit temperature for 
the dominant component. The error in the determination of 
metallicity for HD 98800 is dominated by systematic effects 
a nd is of the order of 01 d ex, the origin of which is discussed 
in lSoderblom et all J2009). This is not a measure of the inter- 
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TABLE 5 

Derived relative radii versus [M/H] 





Tas, 




^Bb 








[M/H] 


(K) 


(K) 


(K) 




^Bb/RAa 


x 2 


-0.6 


4350 


4050 


3500 


0.71 


0.32 


16.5 


-0.4 


4350 


4200 


3500 


0.71 


0.32 


12.3 


-0.3 


4500 


3900 


3500 


0.80 


0.36 


11.4 


-0.2 


4500 


4350 


3500 


0.80 


0.36 


11.3 


-0.1 


4650 


3750 


3500 


0.80 


0.36 


11.8 


0.0 


4650 


3750 


3500 


0.80 


0.36 


12.5 



nal error in this analysis, which is ~ 0.03 dex (from Tabled 
and negligible compared to the systematics. 

3.6. Flux Ratios and photometry matching 

We find that the line ratios for the three stars change from 
night to night, with the lines for Ba becoming relatively 
stronger from the first to the second night and the lines from 
Aa becoming relatively weaker. We infer (Section 13.31 ) that 
the stars may not have been well placed in the center of the 
spectrograph slit during the observations and hence there is 
some uncertainty in the determination of the true flux ratios 
(at ~ 6090 A) of the stars and hence in their radii. However, 
this does not affect the ratio Bb/Ba. We find the mean and 
standard deviations of this ratio for the lowest 40 \ 2 values 
for the three observing nights to be 0.376 ± 0.007 (Ax 2 = 
0.12), 0.430 ± 0.007 (A x 2 = 0.23) and 0.435 ± 0.010 
(Ax 2 = 0.20). This corresponds to a SNR-weighted Bb/Ba 
average flux ratio of 0.416 ± 0.005. The result does not de- 
pend strongly on the number of low-% 2 solutions averaged. 

HST WFPC2 photo metry for the system carried out by 
ISoderblom et alj dl998l) yielded flux ratios for the components 
with Aa:Ba:Bb::0.62:0.3 16:0.064 in the V band, with Bb/Ba 
of 0.207. We carried out a set of runs by fixing the flux ratios 
of the three stars at these values and fitting the spectrum on the 
first night for temperatures and metallicities. The results are 
shown in Table [5] Using flux ratios based on the photometry 
alone, we find that the best fits favor sub-solar metallicity for 
the system. A value of [M/H] ~ -0.2 to —0.3 is indicated 
from these fits. 

4. CONCLUSIONS 

We have fitted high resolution spectra of the HD 98800 sys- 
tem to determine an ov erall metallicity of — 0.20 ± 0.10. This 
confirms indications bv lBoden etal.1 d2005l) that HD 98800 is 
a system with lower metal abundance than the Sun. Coupled 
with age estimates of HD 98800 placing it between 8 and 20 
Myrs, this would imply that there is ongoing star formation 
in the solar neighborhood forming relatively metal-poor sys- 
tems even in the pr esent epoch. This result is consistent with 
ISantos et all (120081) . who find -0.13 < [ Fe/H] < -0 . 05 for 
six nearby star-forming regions and with iJames et alj (2006) 
who find -0.17 < [Fe/H] < -0.02 for three southern star- 
forming regions. 

We have also estimated the effective temperatures of the 
visible components to be T Aa = 4530 ± 80 K (stat), T Ba = 
4200 ± 90 K (stat), T Bb - 3650 K and the flux ratio for Bb/Ba 
to be 0.416 ± 0.005 at 6090 A. Finally, we have determined a 
new orbital period of 314.89 ± 0.03 days for the B-system. 



We note that it is possible for the apparent metallicity de- 
rived from a stellar spectrum to be lowered in the presence 
of very high activity levels, due to filling in of line cores. 
The HD 98800 system is very young, but the activity seen 
in its components is modest. For example, none of the com- 
ponents exhibits Ha emission strong enough to rise above 
the continuum. The activity seen in HD 98800 is also less 
than in Pleiades m embers of similar color, and, as we show in 
SoderblometaD <|2009), we can measure metallicity in those 
Pleiads without any apparent effect from activity. In partic- 
ular, the abundances show no trend with color in going from 
F dwarfs into G stars in the Pleiades ( Soder blom et al.l 2009). 
despite a systematic rise in activity levels. We conclude that 
our determination of metallicity in HD 98800 is unlikely to be 
affected significantly by stellar activity. 

One potential source of systematic error which could cause 
our metallicity estimate to be biased towards low metallicity 
is our assumption that the faint fourth companion Ab con- 
tributes no light to our spectra. If, instead, Ab had signifi- 
cant flux in our wavelength range (and particularly if it were 
a rapid rotator), it would raise the continuum level and hence 
reduce the measured equivalent widths of all lines - which our 
modeling would interpret as indicative of a lower metallicity. 
Because Ab has been detected now in high-resolution imag- 
ing (B05), we can estimate the effect of Ab on our spectra. 
Boden (private communications, 2008) estimates that at the 
K band, the Ab/Aa contrast ratio is 15%. From the measured 
near-IR photometry for HD98800 and a nominal distance of 
40 pc (B05), we therefore estimate absol ute magnitudes of 
K(Aa) = 3.4 and K(Ab) = 5.5. Using the lSiess el tail (120001) 
10 Myr isochrones, these absolute K magnitudes correspond 
to spectral types about K6 and M6. The Siess isochrones also 
provide estimated broadband R — K colors, with V — R ~ 2.3 
for K6 and V — R > 5.7 for M6, yielding R-band magnitudes 
of i?(Aa) = 5.7 and R(Ab) = 11.3. This indicates that the 
flux ratio of Ab/Aa in the wavelength range from which we 
derive metallicities is less than 1%; the influence of Ab on 
the combined spectrum including Ba and Bb is even less. We 
therefore conclude that excluding Ab from our analysis has a 
negligible (less than 1 %) influence on our derived metallicity. 

The stellar parameters we determine could be checked for 
consistency with pre-main sequence stellar models by com- 
puting evolutionary tracks with the above parameters and 
checking whether the stars truly lie on a single isochrone. 
Fully resolved separate spectra of the components of HD 
98800 would provide more accurate parameters, which should 
be feasible with adaptive optics systems given the separations 
of the stars. Astrometry using adaptive optics would also im- 
prove estimates of dynamical masses for the components, en- 
abling tighter constraints to be placed on the age of the sys- 
tem. HD 98800 will continue to be an important reference 
system for understanding pre-main sequence evolution. 



We thank G. Torres of the Center for Astrophysics for pro- 
viding information on the orbits of the components of HD 
98800 and A. F. Boden of the California Institute of Technol- 
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auspices of the Summer Student Program at the Space Tele- 
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